Alternating Direction Implicit Method for Two-dimensional Fokker-planck Equation of Dense Spherical Stellar Systems
نویسندگان
چکیده
The Fokker-Planck (FP) model is one of the commonly used methods for studies of the dynamical evolution of dense spherical stellar systems such as globular clusters and galactic nuclei. The FP model is numerically stable in most cases, but we find that it encounters numerical difficulties rather often when the effects of tidal shocks are included in two-dimensional (energy and angular momentum space) version of the FP model or when the initial condition is extreme (e.g., a very large cluster mass and a small cluster radius). To avoid such a problem, we have developed a new integration scheme for a two-dimensional FP equation by adopting an Alternating Direction Implicit (ADI) method given in the Douglas-Rachford split form. We find that our ADI method reduces the computing time by a factor of ∼ 2 compared to the fully implicit method, and resolves problems of numerical instability. Subject headings: stellar dynamics — methods:numerical — globular clusters:general
منابع مشابه
Solving Two-Dimensional Fuzzy Partial Dierential Equation by the Alternating Direction Implicit Method
متن کامل
An Alternating Direction Implicit Method for Modeling of Fluid Flow
This research includes of the numerical modeling of fluids in two-dimensional cavity. The cavity flow is an important theoretical problem. In this research, modeling was carried out based on an alternating direction implicit via Vorticity-Stream function formulation. It evaluates different Reynolds numbers and grid sizes. Therefore, for the flow field analysis and prove of the ability of the sc...
متن کاملA Heterogeneous Alternating-direction Method for a Micro-macro Dilute Polymeric Fluid Model
We examine a heterogeneous alternating-direction method for the approximate solution of the FENE Fokker–Planck equation from polymer fluid dynamics and we use this method to solve a coupled (macro-micro) Navier–Stokes–Fokker–Planck system for dilute polymeric fluids. In this context the Fokker–Planck equation is posed on a high-dimensional domain and is therefore challenging from a computationa...
متن کاملNumerical Studies and Simulation of the Lower Hybrid Waves Current Drive by using Fokker – Planck Equation in NSST and HT-7 Tokamaks
Recent experiments on the spherical tokamak have discovered the conditions to create a powerful plasma and ensure easy shaping and amplification of stability, high bootstrap current and confinement energy. The spherical tours (ST) fusion energy development path is complementary to the tokamak burning plasma experiment such as NSTX and higher toroidal beta regimes and improves the design of a po...
متن کاملAnalysis and Implementation of Numerical Methods for Simulating Dilute Polymeric Fluids
Analysis and Implementation of Numerical Methods for Simulating Dilute Polymeric Fluids David Knezevic Doctor of Philosophy Balliol College Michaelmas Term 2008 In this thesis we develop, analyse and implement a number of numerical methods for simulating dilute polymeric fluids. We use a well-known model in which the polymeric fluid is represented by a suspension of dumbbells in a Newtonian sol...
متن کامل